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We study the evolution of a V-type three-level system, whose two resonances are coherently ex¬ 
cited and coupled by two ultrashort laser pump and probe pulses, separated by a varying time delay. 
We relate the quantum dynamics of the excited multi-level system to the absorption spectrum of 
the transmitted probe pulse. In particular, by analyzing the quantum evolution of the system, we 
interpret how atomic phases are differently encoded in the time-delay-dependent spectral absorption 
profiles when the pump pulse either precedes or follows the probe pulse. We experimentally apply 
this scheme to atomic Rb, whose fine-structure-split 5s 2 Si/2 —► 5p 2 Pi/2 and 5s 2 Si/2 —> f>p 2 P?,/2 
transitions are driven by the combined action of a pump pulse of variable intensity and a delayed 
probe pulse. The provided understanding of the relationship between quantum phases and ab¬ 
sorption spectra represents an important step towards full time-dependent phase reconstruction 
(quantum holography) of bound-state wave-packets in strong-field light-matter interactions with 
atoms, molecules and solids. 

PACS numbers: 32.80.Qk, 32.80.Wr, 42.65.Re 


A long-standing dream in the physics of light-matter 
interaction is the observation of the quantum dynamics 
of electrons bound in atoms and molecules [1-3]. Tra¬ 
ditional spectroscopy methods enable one to access the 
evolution of the amplitudes, related to the probability 
that a quantum system is in a given state. Much more 
demanding is the extraction of the associated quantum 
phases. 

With the development of femtosecond and attosecond 
laser pulses, transient-absorption spectroscopy (TAS) has 
emerged as a powerful method to gain amplitude and 
phase information on the evolution of a quantum sys¬ 
tem [4-9]. In such pump-probe geometry [Fig. 1(a)], a 
short pump pulse of tunable intensity excites a quantum 
system, whose ensuing dynamical evolution is measured 
by a weak probe pulse. The absorption spectrum of the 
transmitted field results from the interference of the in¬ 
coming probe pulse with the electric field emitted by the 
system [10]. Since this interference pattern depends on 
the state of the system at the arrival of the probe pulse, 
its quantum-mechanical evolution can be time resolved 
by modifying the inter-pulse delay. Furthermore, in con¬ 
trast to most strong-field studies relying on the ioniza¬ 
tion of the system [11, 12], TAS enables direct access to 
bound-state dynamics. 

Very recently, experiments started to investigate the 
absorption spectrum of a probe pulse which excites a co¬ 
herent superposition of quantum states and is followed 
by a strong laser field [7, 13, 14]. For nonautoionizing 
bound states, in the absence of this subsequent pulse, 
well known symmetric Lorentzian line shapes are ob¬ 
served. However, an additional strong pulse can be em¬ 
ployed to nonlinearly drive the state generated by the 
probe field, thus manipulating the absorption lines into, 



Figure 1. (Color online) (a) Transient-absorption- 
spectroscopy setup in which, upon interaction of an atomic 
ensemble with two pulses separated by a time delay r, the ab¬ 
sorption spectrum of the transmitted probe pulse is measured, 
(b) V-type three-level scheme used to model the ensemble of 
Rb atoms interacting with the laser pulses. 


e.g., asymmetric Fano-like shapes [15-17]. Although this 
method opens interesting prospects for line-shape quan¬ 
tum control 1 18, 19], it still remains an open question 
how phase information is encoded in the time-delay- 
dependent strong-field absorption spectrum. 

Here, we investigate theoretically and experimentally 
the phase information which can be extracted from 
transient-absorption spectra for both above-mentioned 
scenarios, i.e., when the dynamics are either triggered 
or modified by an intense pump pulse. We use the V- 
type three-level scheme depicted in Fig. 1(b) to describe 
an ensemble of Rb atoms, with fine-structure-split ex¬ 
cited states. From the ground state, the single-electron 
transitions 5 s 2 S'i/2 -A 5p 2 P\/2 (|1) —> |2) at 794.76nm) 
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and 5 s 2 S'i /2 —> 5p 2 P 3 / 2 (|1) —> 13) at 780.03nm) are ex¬ 
cited by fs-duration pump/probe optical pulses centered 
at 780 nm and separated by a time delay t. We focus 
on the oscillations in r exhibited by the absorption spec¬ 
trum [20] and show how these oscillating features differ¬ 
ently encode the information on the intensity-dependent 
atomic phases for positive and negative time delays. 

The atomic system is described by the density matrix 
p(t,r) [21], of elements pij(t,r), i, j £ {1, 2, 3}. An 
off-diagonal element pij(t,r), i ^ j, represents the co¬ 
herence between states |«) and \j). while a diagonal el¬ 
ement pn(t,r) quantifies the population of level |i) at 
energy u,, with transition energies = u>i — uJj given 
by u> 2 i = 1.56 eV and w 3 i = 1.59 eV. Here, atomic 
units are used unless otherwise stated. The probe pulse 
is centered at t = t and is modeled by the classical 
field £ pr (t,r) = £ pT , 0 f{t - r) cos[w L (t - r) + </ pl \e z , 
with the unit polarization vector e z along the z di¬ 
rection, the laser frequency ul = 1.59 eV, the initial 
phase <j > pr = 0, and the peak field strength £ pr ,o = 
yj8nal pr related to the peak intensity / pr and the fine- 
structure constant a. The envelope function is mod¬ 
eled by f(t) = cos 2 (nt/T) R{t/T), where R(x) = 6(x ± 
1/2) — 6(x — 1/2) is defined in terms of the Heaviside 
step function 6(x). Consequently, /(f) is nonvanishing in 
an interval of duration T = 7tTfwhm / [2 arccos ( \J\ /2)], 
with Tfwhm = 30 fs being the full width at half maxi¬ 
mum of / 2 (f). Similarly, the pump pulse is modeled by 
£ pu (t) = £ pu ,o fit) cos(wi/ + /> pu )e 2 , with initial phase 
</>pu = 0, peak field strength £ pu ,o> and peak intensity J pu . 
The pump pulse is centered at t = 0, such that it precedes 
the probe pulse for r > 0 (pump-probe) and follows it for 
r < 0 (probe-pump). Both linearly polarized pulses cou¬ 
ple the ground state |1) = 5s 2 Sl/2, with magnetic quan¬ 
tum number M = ±1/2, to the two closely lying states 
|2) = 5p 2 Pi/2 and |3) = 5p 2 P 3 /2, also with M = ±1/2 
[22], For the pulse intensities used, we consider only 
electric-dipole-(£l-)allowed transitions with AM = 0, 
for which the dipole-moment matrix elements are d\ k = 
d±k e z > k £ {2, 3}, with d 13 = di2V / 2 [20, 23]. This allows 
us to introduce the two complex, time-dependent Rabi 
frequencies f2 Rk (t,r) = d lk [£ pr ,o f(t - t) e~ luJLT e 1 ^ ± 
£ P u,o /(i)e I0pu ], such that d lk • [5 pu (f) ± £ pl (t,r)\ = 
[f2 Rk (t, r)e i “ L * ± f2^ k (t, r)e" iWLt ]/2. 

The atomic system is initially in its ground state, 
i.e., pij : o = SnSji, with the Kronecker symbol 
Su'. The dynamical evolution of the system results 
from the Liouville-von Neumann equation dp/ dt = 
-i[H 0 + flint (i, t), p(t, r)[ - i{P, p(t, r)}, with [A, B} = 
AB — BA and {A,B} = AB ± BA [21]. Here, 
flo = y/?_ 1 oji&u is the electronic-structure Hamil¬ 
tonian, with ladder operators dtj = |i)(/|, while 

flint(f,r) = - 1=2 ai k e lulIjt + Ji.c. is the 

fll-light-matter-interaction Hamiltonian in the rotating- 
wave approximation. The diagonal matrix r has ele¬ 


ments Pjj = Sijjj, with 71 = 0 and 72 = 73 = 1/(500fs). 
These decay rates are much larger than the spontaneous 
decay rates of the system and effectively account for 
broadening effects at work in the experiment. 

The measured optical-density (OD) absorption spec¬ 
trum is given by ^(w, r) = - log[5 pr!OUt (a;, r)/5 pr ,in(w)], 
where 5 pr> i n (a;) and 5 priOU t(w, r) are the spectral intensi¬ 
ties of the incoming and transmitted probe pulse, respec¬ 
tively. The absorption spectrum is related to the single¬ 
particle dipole response of the system, i.e., the evolution 
of the atomic coherences pvi{t. t) and pi 3 (t,r), via [ 10 ] 


(w, t) oc — Im 



u k \\di k |pifc(t, r)e 1Uj(t r) dt 

(1) 


with Fourier transforms centered at the arrival time 
of the probe pulse. In the presence of a single 
pulse, the system, excited at t = to into the state 
p(to), decays with time-dependent coherences pi k (t) = 
Pifc(io) e 1Wfcl ( t-t °) giving rise [Eq. (1)] to sta¬ 

tionary lines centered on the transition energies 1 x 21 and 
w 3 i, with linewidths given by the decay rates 72 and 7 3 , 
respectively. However, if two pulses are used, the coher¬ 
ences pife(t, r) become explicitly time-delay-dependent, 
implying a modulation in r of intensity and shape of the 
spectral lines. Although rapid oscillations in r character¬ 
ize the absorption spectrum in Eq. (1), these fast oscil¬ 
lating features are averaged out for the noncollinear ex¬ 
perimental geometry in Fig. 1(a). We include this effect 
in our model via convolution of the spectrum in Eq. (1) 
with a normalized Gaussian function G(r, At) of width 
Ar = 5 x 27t/u;l, thus obtaining 


/ OC 

G(t—t\ At) Sf l(w, t') dr'. 

-OO 

( 2 ) 

In order to investigate the intensity-dependent phase 
information which can be extracted from J^(oj,t), in 
Fig. 2 we compare numerical calculations of Eq. (2) for 
a probe intensity of 9.2 x 10 8 W/cm 2 and different pump 
intensities. As shown in Fig. 2(a) for a peak intensity of 
1.0 x 10 10 W/cm 2 , the spectrum consists of two lines cen¬ 
tered on the fll-allowed-transition energies ui 21 and w 3 i. 
Furthermore, for both positive and negative time delays, 
amplitude and shape of the two absorption lines oscil¬ 
late in |r|, with a periodicity determined by the beating 
frequency w 3 2 = 27r/(140fs) and a discontinuity across 
r = 0, where the pulses overlap. The black lines in 
Fig. 2(a) highlight this behavior, by displaying ^(1x21, t) 
and 5^ (w 3 i, r) evaluated at the atomic transition ener¬ 
gies. Absorption strengths are shown in arbitrary units, 
as they depend on the gas pressure in the cell, which is 
not included in the simulation. 

Figure 2 (b) displays analogous features for a pump in¬ 
tensity of 2.8 x 10 lo W/cm 2 . Here, however, different 
phases of the time-delay-dependent oscillations can be 
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Figure 2. (Color online) Theoretical absorption spectra of a 
delayed probe pulse of intensity 9.2 x 10 s W/cm 2 interacting 
with an ensemble of Rb atoms, modeled via the three-level 
system in Fig. 1(b) and driven at t = 0 by a pump pulse of 
intensity (a) 1.0 x 10 10 W/cm 2 and (b) 2.8 x 10 10 W/cm 2 . 
Panels (c) and (d) display, for negative and positive time 
delays, respectively, absorption spectra for a pump pulse of 
1.0 x 10 10 W/cm 2 and artificial phases (j >2 and (p 3 [Eq. (3)], 
with (c) <f >2 = 4>3 = —2.26 for r < 0, and (d) (p 2 = 0, (f >3 = 0.67 
for r > 0. In each panel, the top (bottom) black lines repre¬ 
sent the absorption spectra evaluated at the transition energy 
CJ 31 ( 0 J 21 ) in arbitrary units. All black lines are on the same 
scale, with the 0 aligned on the corresponding transition en¬ 
ergy. 


recognized. This is apparent from the time delays at 
which o 5 ^(w 2 i,r) and ,r) exhibit, e.g., their min¬ 

ima, as highlighted by the red, dashed lines. For r > 0, 
these minima are aligned both in Fig. 2(a) and 2(b), al¬ 
though a common phase shift can be distinguished. In 
contrast, for r < 0, the minima are aligned in Fig. 2(a), 
while in Fig. 2(b) they appear shifted towards opposite 
directions as a consequence of the larger pump intensity. 
The different behavior displayed at positive and negative 
delays reflects the different excitation sequence under¬ 
gone by the system. While at r > 0 the pump pulse 
generates the initial coherences, for r < 0 it strongly 
modifies the coherent superposition of states generated 
by the probe pulse. 


In Figs. 2(c) and 2(d) we show that the phase of the 
oscillations in r encodes the key information about the 
intensity-dependent phases of the atomic system. In both 
panels, a pump intensity of 1.0 x 10 10 W/cm 2 is used as 
in Fig. 2(a). Here, however, we modify the phases of 
the atomic system by artificially transforming the density 
matrix p(T/2,r)—describing the state of the system at 
t = T/2 immediately after the interaction with the pump 
pulse—into 

p'(T/2,r)=£p(T/2,r)£l (3) 

The usual equations of motion are employed to calculate 
the subsequent dynamics of the system. In Eq. (3), £ 
is a diagonal matrix of elements r I> t j = e _1 ^ i Sij , with 
cf> i = 0 , while <f>2 and cp 3 represent shifts in the phase of 
levels 12) and |3), respectively, such that p'^iT/2, r) = 
pij(T/ 2 ,T)e-W*-M. 

These artificial shifts (pk effectively model phase 
changes undergone by the system at varying pump in¬ 
tensities, to which the time-delay-dependent absorption 
spectra are sensitive. For instance, in Fig. 2(c), each 
(pk is set equal to the difference between the phase 
changes caused by a pump intensity of 2.8 x 10 lo W/cm 2 
[Fig. 2(b)] and of 1.0 x 10 10 W/cm 2 [Fig. 2(a)], which can 
be extracted from the corresponding dynamics of the den¬ 
sity matrix. The resulting absorption spectrum features 
time-delay-dependent oscillations with the same ampli¬ 
tudes as in Fig. 2(a), but whose phases are now equal to 
those displayed in Fig. 2(b). An analogous procedure is 
applied in Fig. 2(d) for positive time delays, with new 
phases (pk- 

In order to test the relation between intensity- 
dependent atomic phases and time-delay-dependent ab¬ 
sorption spectra, in Fig. 3 we show results from the ex¬ 
perimental measurement of the OD absorption spectrum 
of a probe pulse transmitted through an ensemble of Rb 
atoms. A Ti:Sapphire laser is utilized, with a repetition 
rate of 4 kHz, pulse energy of 0.7 mJ, FWHM of 30 fs, and 
central wavelength of 780 nm. The laser beam is split by 
a spatial mask with two irises into two beams, which are 
used as probe and pump beams in Fig. 1(a). A delay 
line is included into the pump beam path by means of a 
piezo translation stage at close-to-normal incidence, al¬ 
lowing one to scan the time delay between the two pulses. 
Rb vapor is supplied in a cell with 20 mm length, which 
is heated to 160 °C by a home-made heating and con¬ 
trol system. Using the known vapor pressure of Rb, 
the atomic density is estimated to be 2.8 x 10 14 cm -3 . 
The pump and probe pulses, with beam sizes of 3.5 mm 
and 3.0 mm, respectively, are focused by a concave mirror 
with a focusing length of 500 mm. After passing through 
the sample, the probe beam is picked up by an identical 
concave mirror and coupled into a fiber-pigtailed spec¬ 
trometer (Ocean Optics, USB2000+). The probe pulse 
has an intensity of 9.2 x 10 8 W/cm 2 , and Figs. 3(a) and 
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Figure 3. (Color online) Experimental absorption spectra of 
a delayed probe pulse of intensity I pl = 9.2 x 10 8 W/cm 2 , 
transmitted through an ensemble of Rb atoms driven at t = 0 
by a pump pulse of intensity (a) I pu = 1.0 x 10 10 W/cm 2 and 
(b) I pu = 2.8 x 10 10 W/cm 2 . The black on-resonance lineouts 
have the same meaning as in Fig. 2. 


3(b) show experimental results for a pump intensity of 
1.0 x 10 10 W/cm 2 and 2.8 x 10 10 W/cm 2 , respectively. 
The absorption lines at each resonance are broadened to 
a measurable width as a result of Doppler and collision- 
induced pressure broadening, as well as additional broad¬ 
ening in the nanosecond pedestal of the femtosecond 
pulses. For low pump intensities [Fig. 3(a)], at r < 0, the 
small amplitude of the oscillations in r renders the asso¬ 
ciated phases not straightforwardly extractable. Apart 
from that, the experimental results are in good agree¬ 
ment with the theoretical predictions from Fig. 2, show¬ 
ing analogous intensity-dependent phase shifts in the os¬ 
cillating features of the absorption spectra. 

Finally, we focus on the dependence of phase effects 
upon the intensity of the pump pulse. Figures 4(a) 
and 4(b) display the experimental absorption spectra 
o 5 ^(w 2 i,t) and ^(w 3 i,r) at the transition energies ui 2 i 
and W 31 , respectively, for varying pump intensities. These 
measurement results quantify intensity-dependent shifts 
in the phase of the time-delay-dependent spectra, which 
are in good agreement with the theoretical predictions 
in Figs. 4(c) and 4(d) and confirm the observations from 
Figs. 2 and 3. 

The experimentally observable phase properties dis¬ 
played by J^(w 21 ,t) and ^(w 3 i,r) reflect the different 
physical mechanisms acting at negative and positive time 
delays. This is further analyzed by implementing the 
artificial-phase method introduced in Eq. (3) to repro¬ 
duce intensity-dependent phase effects. For r < 0, the 



Phase (rad) Time delay (fs) Phase (rad) 


Figure 4. (Color online) [(a)-(b)] Experimental and [(c)-(d)] 
theoretical absorption spectra at fixed frequencies [ 0 J 21 = 
1.56eV in (a) and (c), LO 31 = 1.59 eV in (b) and (d)] of a de¬ 
layed probe pulse of intensity Ipr = 9.2 x 10 8 W/cm 2 , trans¬ 
mitted through an ensemble of Rb atoms driven at t = 0 
by a pump pulse of varying intensity. The dashed, magenta 
lines display the time delays (as a function of pump inten¬ 
sity) associated with minima in the absorption spectra. Pan¬ 
els (e), (f) and (g) show the artificial phases \ 4 >2 in (e), <j >3 
in (f) and (g), see also Eq. (3)] which need be added at 
7 pu = 1.0 x 10 10 W/cm 2 to reproduce the intensity-dependent 
phases of the time-delay-dependent oscillations at [(e)-(f)] 
negative and (g) positive time delays. 


pump pulse strongly modifies the excited state generated 
at t — t by the probe pulse. The intensity-dependent 
phases of the oscillations in r featured in Figs. 4(c) and 
4(d) enable one to quantify the shifts in the phase of 
the atomic coherences ensuing from the interaction with 
strong pump pulses of varying intensity, as captured by 
<h(Ipu) an d </> 3 (/pu) in Figs. 4(e) and 4(f), respectively. 
In contrast, for r > 0, the pump pulse is responsible for 
the initial excitation of the system. The parallel shift un¬ 
dergone by c 5 ^(w 2 i, t) and ^( 0231 , r) at increasing pump 
intensities [Figs. 4(c) and 4(d)] is therefore a consequence 
of the dependence of the initial coherences—and their rel¬ 
ative phase—upon the intensity of the pump pulse. Ar¬ 
tificial phases can be used to effectively reproduce this 
behavior. In this case, however, the phases of the re¬ 
sulting spectra ^(uj 2 i,t) and = 5 ^(w 3 i,t) are identically 
affected only by the phase difference </>3 (7 pu ) — §2 (/ pu ), 
which is shown in Fig. 4(g) for (f >2 (Ipu) = 0. 

In conclusion, we have used a V-type three-level sys¬ 
tem to investigate strong-field-induced phase effects in 
transient-absorption spectra for both positive and neg¬ 
ative time delays and for pump pulses of variable in¬ 
tensity. Furthermore, we have developed an artificial- 
phase method to interpret time-delay-dependent absorp- 
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tion spectra and relate them to the phase evolution 
and intensity dependence of atomic quantum coher¬ 
ences. Although this method was here introduced to an¬ 
alyze intensity-dependent phase effects, the level-specific 
phases (f>k could also be employed to model additional 
interactions of the atomic system with external fields se¬ 
lectively modifying the two excited states. Understand¬ 
ing the phase-sensitive link between quantum dynamics 
and transient-absorption spectra represents an essential 
step towards the application of TAS to the extraction 
of quantum phases in bound-state atomic systems inter¬ 
acting with strong-field pulses and to sensitive tests of 
quantum-dynamics theory. 
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